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We discuss a new leading-order parameterization of generalized parton distributions of the pro- 

; 

, ton, which is based on the idea of duality. In its minimal version, the parameterization is defined 

in ' 

by the usual quark singlet parton distributions and the form factors of the energy-momentum 



tensor. We demonstrate that our parameterization describes very well the absolute value, the 
Q 2 -dependence and the VP-dependence of the HERA data on the total DVCS cross section and 
contains no free parameters in that kinematics. The parameterization suits especially well the 
low-x^j region, which allows us to advocate it as a better alternative to the frequently used double 
distribution parameterization of the GPDs. 
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I. INTRODUCTION 



Generalized parton distributions (GPDs) have become a standard QCD tool for ana- 
lyzing and parameterizing the non-perturbative parton structure of hadronic targets, for 
reviews see 0, 0, S, , S, la| . In general, GPDs are more general and complex objects than 
structure functions and form factors. In addition, experimentally measured observables do 
not access the GPDs themselves but only their convolution with the hard scattering coef- 
ficients. Therefore, the experimental determination of the GPDs is an extremely difficult 
task. Hence, when dealing with the GPDs, one invariably uses models, the known limiting 
behavior and general properties of the GPDs and the physical intuition. 

The GPDs have been modeled using virtually all known models of the nucleon structure: 
bag models M, the chiral quark-soliton model lg, light-front models M, llCI, constituent 

[ n n 

quark models 11 , Bethe-Salpcter approach 12 , a NJL model 13 . In addition, a double 

nfi n 

distribution model of the GPDs [1J, [15|| and modeling by perturbative diagrams [16| were 
suggested. 

The factorization theorem for deeply virtual Compton scattering (DVCS) gives a 
practical possibility to measure the GPDs by studying various processes involving the GPDs: 
DVCS, exclusive electroproduction of vector mesons, wide angle Compton scattering r| 
lig| | . exclusive pp — > 77 annihilation j^j], the pp — > 771- process J21], 7*7 — > mr near 



threshold 



22j. However, in order to accommodate such a potentially large number of data, 



parameterizations of the GPDs should be sufficiently flexible and versatile. In particular, 
they should allow for the connection of the DVCS with the pp — > 7 7 process. 

The commonly used double distribution parameterization of the GPDs ^| is one 
example of the model of the GPDs which could be used to connect different physical chan- 
nels [23]. However, the parameterization of the GPDs based on the double distribution has 
a number of problems from the phenomeno logical point of view. First, in order to have the 
full form of polynomiality, the so-called D-term |24| has to added by hand. Second, the 
model dramatically overestimates the low-x HERA data on the total DVCS cross section 
because it involves proton parton distributions at extremely small and unmeasured values 
of Bjorken x Third, the model does not allow for an intuitive physical motivation and 
interpretation, see 2(| for a discussion of the physics of GPDs. 

In this paper, we offer a new model for the GPDs, which was in a general form intro- 
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duced in |2j|. Unlike the models of the GPDs mentioned above, the present model has a 
simple phy sical interpretation and direct correspondence to the mechanical properties of the 
target [28[ . The suggested parameterization of the GPDs can be analytically continued in t 
to the physical region of the pp — > 7 7 reaction and also allows for flexible modeling of the 
t-dependence of the GPDs. 

The considered parameterization of the GPDs is called dual because the GPDs are pre- 
sented as an infinite series of t-channel exchanges, which reminds of the ideas of duality in 
hadron-hadron scattering. 

In this work, we formulate the minimal version of the dual parameterization and deter- 
mine the free parameters of the model. Using the resulting dual parameterization of the 
GPDs, we successfully describe the HERA data on the DVCS cross section 

0, Q Q. 

We explain that our parameterization suits the low-x Bj kinematics especially well because 
the quark singlet parton distributions are never probed at the unmeasurably low values of 
Bjorken x and because the final expression for the DVCS amplitude is numerically stable. 
Thus, the dual parameterization of the GPDs gives an opportunity to have a physically 
intuitive model of the GPDs, which agrees with the DVCS experiments and which can serve 
as an alternative to the popular double distribution model. 



II. THE DUAL PARAMETERIZATION OF GPDS 

The dual representation of the GPDs was first introduced for the pion GPDs in (3^] . The 
essence of that derivation is presented below. The starting point was the decomposition 
of the two-pion distribution amplitude (27rDA) in terms of the eigenf unctions of the ERBL 
evolution equation (Gegenbauer polynomials C^ 2 ), the partial waves of produced pions 
(Legendge polynomials Pi) and generalized form factors B n \. The moments of the 27rDA, 
being the matrix elements of certain local operators, could be related by crossing to the 
moments of the pion GPDs. Then, the pion GPDs could be formally reconstructed using 
the explicit form of their moments. 



in 



Based on the result of 32], the dual representation for the proton GPDs was suggested 
In this paper, we will consider only the singlet (C-even) combination of the GPDs 
H, which give the dominant contribution to the total DVCS cross section at high energies 
and small t. We will work in the leading order approximation and, hence, we will consider 
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only quark GPDs. 

The dual representation of the singlet GPD H l of the quark flavor i is [2 



h\x,u^ 2 ) = E E But^ 2 )e (1 - (1 - pi cT (?) Pi (±) , (i) 

n=l, odd i=0, even \ 



where x, £ and t are the usual GPD variables. The series (JIJ is divergent at fixed x and 
£, and, hence, it should be understood as a formal series. In particular, it is incorrect to 
evaluate the series term by term. As a result, the GPD H % of Eq. (0) has a support over the 
entire — 1 < x < 1 region, regardless that each term of the series is non- vanishing only for 
— £, < x < £. The formal represent at ionO can be equivalently rewritten as a converging 
series using the technique developed in 27| . 

The derivation of Eq. (0) used the idea of duality of hadronic physics, when the scattering 
amplitude in the s-channel is represented as an infinite series of the t-channel exchanges. 
This explains the name "the dual representation" for the representation of Eq. 

As a double series, Eq. Q is inconvenient for phenomenological applications. For the 
evaluation of the LO DVCS amplitude, it is useful to introduce the functions Qk(x, t) whose 
Mellin moments generate the B l nl form factors 27] 

I^ +1 _ fe (t,/i 2 ) = f dxx n Ql(x,t^ 2 ), (2) 
Jo 

where k is even. A remarkable property of the dual representation is that the /^-evolution 
of the functions Q\ is given by the usual leading order (LO) DGLAP evolution. 

Since the B l nl form factors are related to the moments of H l , the Qk functions are also 
constrained by these moments. In particular, from 



[ l dx x H\x, e, t) = M*(t) + U 2 d l (t) = ^ 
J-i 5 5 



B 12 (t)-^(B 12 (t)-2B 10 (t))e 



(3) 



it follows that 



/ dxxQ {x,t,fi 2 ) = -M*(t,/i 2 ), 
Jo o 

f dx x Q^x, t, fi 2 ) = A Ml(t, fi 2 ) + d l (t, fi 2 ) , (4) 



where M % 2 at t = is the proton light-cone momentum fraction carried by the quarks; d l (t) 
is the first moment of the quark D-term. 
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In addition, the B nn+ i form factors at the zero momentum transfer are fixed by the Mellin 
moments of the quark singlet parton distribution functions (PDFs). In particular, 



- J dx (V (x, /i 2 ) + cf(x, /i 2 )) = Sio(O) = j dx Ql(x, 0, /j 2 ) . (5) 
The Qq functions at t — are completely fixed in terms of the forward proton PDFs 



QHx, 0, /i 2 ) = q\x, /i 2 ) + ^(x, /i 2 ) ~ X -\ X % Uiz, /x 2 ) + q\z, /i 2 )) . (6) 

Z Jx Z v ' 



As suggested in 



27| , keeping only the functions Qq and Q 2 constitutes the minimal version 



of the dual parameterization of the GPDs. The functions Q and Q 2 are defined by Eqs. (JHJ) 
and (jH), where M|(£) and <i 4 (t) have a clear physical interpretation since they are the form 
factors of the energy-momentum tensor evaluated between the states representing the given 
target. At t = 0, M^O) is the fraction of the plus-momentum of the nucleon carried by 
the quarks of flavor i; d l (0) characterizes the shear forces experienced by the quarks in the 
target. 

Next we discuss the minimal version of the dual representation in detail. While Q l at 
t = is defined by Eq. (jH}, only the first x- moment of Q\ is constrained. We simply assume 
that Q\ oc Qq and take 

Q^O,/! 2 ) =PQ (x,0, f i 2 ), (7) 
where (3 l are constants. From Eq. we obtain 

/3 ~5M|(0) + 2' (8) 

which gives 

(3 U = -AA } /3 d = -8.9 , (3 s = 0.5. (9) 

In this numerical estimate, we assume that d u = d d w -2 and d s w 0, which agrees with 
the SU(3)-symmetric chiral quark soliton model calculation of the nucleon D-term [33J: 
X)id*(0) ~ —4, and takes into account the SU(3) symmetry breaking in the nucleon PDFs 
(the suppression of the strange quark PDF at the low resolution scale). The momentum 
fractions M % 2 were evaluated at /xq = 0.6 GeV using the LO GRV parton PDFs Q- 

In general, (3 l also depend on /i 2 . However, as will be seen from the general expression 
for the DVCS amplitude, at small values of £ typical for the HERA data on the total DVCS 
cross section, the contribution of the Q 2 function is kinematically suppressed. Therefore, 
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the goodness of the description of the data is not affected by the exact values of (3\ and we 
simply used Eq. © at all /i 2 . 

Until recently, the t-dependence of the DVCS cross section was not measured. One would 
simply assume that the DVCS cross section exponentially depends on t, 

It = 6XP { ~ B |t|} { It 

such that the total DVCS cross section is 

/ ^ 1 fda DYC s(xBj,Q 2 ,t)\ 
<xdvcs(^,Q ) = - ^ Jt . 

The value of the slope parameter B was rather uncertain, 5 < B < 9 GeV -2 . The range of 
the values covers the experimentally measured range of the t-slope of electroproduction of 
light vector mesons at HERA. However, very recently the t-dependence of the total DVCS 
cross section for 0.1 < \t\ < 0.8 GeV 2 and at Q 2 = 8 GeV 2 was measured by the HI 
collaboration at HERA and was fitted by the exponential form of Eq. (jlUJl with the result 
B = 6.02 ± 0.35 ± 0.39 GeV" 2 

In our numerical estimates of the DVCS cross section, we calculate the DVCS amplitude 
at t = and then use Eq. (|llj) in order to find the t-integrated DVCS cross section. In 
general, the slope B should decrease with increasing Q 2 . A particular model for the in- 
dependent slope was suggested in [2^: B{Q 2 ) = 8 (1—0.15 ln(Q 2 /2)) GeV -2 . In our analysis, 
we use the same Q 2 -dependence, 

B{Q 2 ) = 7.6 (l - 0.15 ln(Q 2 /2)) GeV -2 , (12) 

but with a slightly smaller constant 7.6 GeV" 2 , which is chosen such that Eq. (fT2"|) reproduces 
the HI value of the slope at Q 2 = 8 GeV 2 . 

In summary, our parameterization of the GPDs H % is defined by Eqs. (J3J), © and (J7J). 
The t-dependence of the DVCS cross section is given by Eqs. (jlOj) and ()12j) . This is the 
minimal version of the dual representation of the GPDs, which can be readily extended by 
considering more Q\ functions, a more elaborate ^-dependence and by taking into account 
the other GPDs of the proton. The main practical advantage of our representation is that 
the /^-evolution of Q l 2 is given by the usual DGLAP evolution of the singlet PDFs, see 
Eq. ©. 
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III. DESCRIPTION OF LOW-x HERA DATA ON DVCS CROSS SECTION 



In this section, we evaluate the total DVCS cross section using the minimal model of the 



dual representation of the GPDs and compare the results to the HERA data 
The total unpolarized DVCS cross section on the photon level reads, see e. 



. g . y, 



(rixvcs{xBj,Q ) = . r = = = / dt\Amc${Z.,t,Q ) , (13) 

O 4 , /I -L- 4. r 2 m 2 JQ2 J t . 



Q^l + 4x%m%/Qi 

where a e . m . the fine-structure constant; £ = l/2xBj/ (1— in the Bjorken limit; t maa; ~ 

and i m j n ~ —1 GeV -2 ; |^4dvcs| 2 is the squared and spin-averaged DVCS amplitude. 

To the leading order in a s , the DVCS amplitude is expressed in terms of the singlet 
combination of the GPDs H l , 



A 



DVCS 



(£, t ,<n= ? ^/ 1 ^tf(,,t,^(_i + _i ) . d4) 



£ + zO a; + £ 
the DVCS amplitude can be presented 



Using our model for the GPDs and the results of 
in a compact form in terms of the Q % and Q l 2 functions 

fc=0 



r 1 rjrv. " I 

(15) 

Using the exponential ansatz for the t-dependence of the DVCS cross section, the total 
DVCS cross section is expressed in terms of the DVCS amplitude at t — [see Eq. (fTTj) ] 

where *Advcs(£ ; t — 0, Q 2 ) is given by Eq. (JT5Jl evaluated with Qq 2 (x, 0, Q 2 ). 

Our predictions for the Q 2 and W^-dependence of the total DVCS cross section are pre- 
sented in Figs. Q and respectively. For comparison, we also present the HI [2jJ and 
ZEUS (30 1 data. 

Note that the ZEUS data points, which were taken at W = 89 GeV and at Q 2 = 9.6 GeV 2 , 
have been rescaled to the HI values of W = 82 GeV and Q 2 = 8 GeV 2 using the fitted W 
and Q 2 -dependence of the DVCS cross section: ctdvcs ^ W ' 75 and odvcs ^ l/iQ 2 ) 1 ' 54 30] . 

For the proton forward PDFs, which are required to evaluate Q l , we used the LO CTEQ5L 
parameterization 3t1 ]. 
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20 40 60 80 



Q 2 (GeV 2 ) 

FIG. 1: The total DVCS cross section at W = 82 GeV as a function of Q 2 . The predictions of the 
dual parameterization of the GPDs (solid curve) is compared to the HI 

y 

and ZEUS The 

error bars represent the statistical and systematic uncertainties added in quadrature. 



One can see from Fig. [I] that the absolute value and the Q 2 -dependence of the total DVCS 
cross section is described very well. The agreement with the data at the highest values of 
Q 2 would have been worse, if we had used the Q 2 -independent slope B. 

From Fig. |21 one can see that the absolute value and the lU-dependence of the DVCS 
cross section is also reproduced rather well. However, one should note a slight discrepancy 
between the ZEUS and HI data points at lower values of W and large experimental errors 
at the high end of W. 

It is important to emphasize that our predictions for the total DVCS cross section were 
made using the parameterization of the GPD, which contains no free parameters (the role of 
Q\ and (3, see Eq. (J7J), is negligible in the HI and ZEUS kinematics). It is very remarkable 
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• HI 

O ZEUS 



Q 2 =8 GeV 2 




20 40 60 80 100 120 140 
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FIG. 2: The total DVCS cross section at Q 2 -- 

dual parameterization of the GPDs (solid curve) is compared to the HI 
error bars represent the statistical and systematic uncertainties added in quadrature. 



GeV function of W . The predictions of the 

Q and ZEUS 30]. The 



that the agreement with the data is so good. 

In order to understand, at least partially, the success of the dual parameterization of 
GPDs in the description of the low-x HERA data, it is instructive to analyze the DVCS 
amplitude ^4dvcs of Eq. (|T5j) in some detail. Evaluating the imaginary and real parts of 
Eq. (|15|). one obtains |27( 

Im g 2 ) = - E t — , 1 E x k Q k (x, 0, Q 2 ) , 

Re t) = - E e 2 / — E xk Qk(x, 0, Q 2 ) ( -, ~ 2S k0 ) 

i Ja x t^o K Jl + 2x/i + x 2 ' 
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2>?/ ™J2x k Q k (x,0,Q 2 )( 




sfl - 2xj i + x 2 y/l + 2x/£ + x 2 



1 1 



24o), (17) 



where a = (1 — \/l — £ 2 )/£. 

Al low scgj, £ ~ and the integration limit is a « £/2 = £_Bj/4. Thus, the functions 

Qq and are never sampled at x < XBj/^- This is clearly an advantage over the double 
distribution parameterization of GPDs, where the forward parton distributions are sampled 
all the way down to x — 0, which results in the acute sensitivity to the unmeasured, very 
\cw-x behavior of the proton PDFs and leads to a gross overestimate of the data 

In addition, Eqs. (JT7j) are convenient for the numerical implementation since the inte- 
grands do not contain large end-point contributions, as can be explicitly seen by changing 
the integration variables. 

IV. DISCUSSION AND CONCLUSIONS 

We presented and discussed a new leading order parameterization of GPDs introduced 
in (27]]. In its minimal form, the parameterization is defined by the forward singlet quark 
PDFs and the form factors of the energy-momentum tensor, see Eqs. (jlj) and (JSJ). The 
t-dependence of the DVCS cross section was assumed in a simple factorized form with the 
Q 2 -dependent slope, see Eqs. (fTUJl and (fT2~J). 

We showed that our parameterization of the GPDs describes very well the absolute value, 
the Q 2 -dependence and P^-dependence of the HERA data on the total DVCS cross section. 
Moreover, since the data is at low Xsj, our parameterization can be simplified by omitting the 
contribution of the Q l 2 function. This means that we achieved a remarkably good description 
of a large set of the data on DVCS using a parameterization of the GPDs which contains no 
free parameters! 

We discuss that our parameterization suits the low-x^j kinematics especially well because 
the quark singlet PDFs are never probed at the unmeasurably low values of Bjorken x, as is 
the case for the popular double distribution model and because the expression for the 
DVCS amplitude is numerically stable. This allows us to advertise our model as a better 
alternative to the popular double distribution parameterization of the GPDs, at least in the 
low-£ region. 

The parameterization presented in this work can be readily generalized by including more 
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Qi functions, considering the GPDs E, H and E and by using more elaborate models of the 

n n 

t-dependence. This was not necessary in the HI [31] and ZEUS [30] kinematics, but might 
be required for the HERMES and CLAS kinematics. 

Also, the role of next-to-leading order (NLO) corrections and higher twist effects should 
be investigated. In particular, it is important to compare the size of the NLO corrections 
using the dual parameterization with the results of the analysis using the double distribution 
parameterization, where the NLO corrections were found to be large 
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